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In this work we use multiple scattering in conjunction with a genetic algorithm to reliably determine 
the optimized photonic- crystal-based structure able to perform a specific optical task. The genetic 
algorithm operates on a population of candidate structures to produce new candidates with better 
performance in an iterative process. The potential of this approach is illustrated by designing a 
spot size converter that has a very low F-number (F=0.47) and a conversion ratio of 11:1. Also, 
we have designed a coupler device that introduces the light from the optical fiber into a photonic- 
crystal-based waveguide with a coupling efficiency over 87% for a wavelength that can be tuned to 
1.5 fim. 
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A new generation of optical devices is envisaged thanks 
to the properties of photonic crystals (PC's).^ Though 
the recent advances in three-dimensional PC's structures, 
in the last years much attention has been focused on 
systems based on two-dimensional (2D) PC's because of 
their easiness in the fabrication process. Thus, very com- 
pact optical devices and circuits can be designed by in- 
troducing point and/or line defects. In order to use such 
PC circuits in actual applications it is necessary to es- 
tablish a connection with an optical fiber. However, the 
core of the optical fiber is about one order of magnitude 
larger than the PC-based waveguide. Therefore, the de- 
sign of an efficient (low loss) spot size converter is a cru- 
cial goal in the field of PC; its solution will introduce the 
PCs devices in the market place. In this regard, several 
groups^"'' have tackled this problem by using different 
approaches. Most of them proposed tapered waveguide 
structures'"^, or by using reflective structures to focus 
the light into the waveguide.^ A different approach con- 
sists of using the anisotropy of the PC's equifrequency 
surfaces.'' 

This letter introduces a method that is useful in de- 
termining the optimized configuration of a 2D-PC struc- 
ture capable of performing a requested optical task with 
high efficiency . The method is illustrated by finding a 
spot size converter (lens) that has a conversion ratio 11:1. 
In addition, the designed PC structure that involves a 
spot-size converter in connection with a PC-based waveg- 
uide it is presented. The insertion loss predicted for this 
new structure is about 13%, which is of the lowest re- 
ported by numerical simulations based on different cou- 
pling mechanisms.'"^ 

Our method is based on a binary-coded genetic algo- 
rithm (GA), an optimization strategy inspired by Dar- 
winian evolution*. This method has been applied to 
solve a wide variety of problems in different fields like, 
for example, molecular geometry optimization^, material 
design^", and artificial intelligence^^. In the field of op- 
tics, the GA has been employed in the synthesis of Bragg 
gratings that conform to a particular spectrum, "'^', phase 



recovering from a fringe pattern^'^, and in designing ir- 
regular lateral tapering.' 

Although our proposal is general and applicable to 
any dimensionality, here we analyze 2D-PCs for simplic- 
ity. We consider hexagonal lattice of Si dielectric cylin- 
ders {nsi «3.46 at 1.5 /im) in a background of silica 
(f^Si02=l-45). We employ the multiple scattering the- 
ory (MST) to study the diffraction effects of the TM 
modes (in-plane magnetic field) on structures based on 
that symmetry. The MST has been successfully applied 
by some of us in the analysis of metallic and dielectric 
clusters based on 2D PCs^-^^i^. 

As the first step in the process, a structure made of a 
fixed number of cylinders is chosen. The GA is then im- 
plemented with the possibility of removing cylinders, but 
it cannot change neither their positions nor their radius. 
A string of binary digits represents a possible structure 
(individual). Each binary digit (1 means presence and 
absence of a cylinder at a given position) is called a gene. 
Each individual is associated with a value, the fitness, 
which is strictly related to the PC property that we are 
looking for. A crystal with a good-asked for property 
will have a high fitness and vice versa. At start, a cer- 
tain number of individuals (population) is created by the 
GA that randomly determines their genes. The size of 
population P depends on the complexity of the problem, 
which is strictly related to the number of genes. 

After the initiation, the GA selects two individuals at 
random, the one with the highest fitness acts as a parent, 
the second parent is chosen in the same way. Now, an 
uniform-crossover operator randomly exchanges the chro- 
mosomes of the two parents with a probability C and, in 
this way, an offspring is created. This is done until we 
have a new equal sized population of offspring, i.e. a sec- 
ond generation. The selection and crossover operators 
tend to enable the evolutionary process to move toward 
promising regions of the search space. In a step further, 
a mutation operator is introduced to prevent premature 
convergence to local maximum. This operator changes 
the value of the genes at random with a probability /i. A 
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high value in the probabihty of mutation indicates that 
the GA looks for better solutions at larger distances from 
a local maximum. The three operators are repeatedly 
applied to the population, with the constrain that the 
best-fitted individual is always copied into the next gen- 
eration (elitism). As there is no way to know if the global 
maximum is reached, this process is iterated until there 
is no improvement of the fitness value of the best indi- 
vidual in the population. For further details the reader 
is referred to books by HoUand^and Golberg.^^ 

The fitness parameter in our calculations is equal to the 
component of the Poynting vector, ~^ , along the propa- 
gation direction of the incident beam, which is evaluated 
in a selected position or in a set of positions depending on 
the desired result. The Poynting vector is calculated by 
the MST. The external radiation field is here a Gaussian 
beam in order to demonstrate the spot-size conversion. 
The beam is formed as a sum of plane waves weighted 
with a function dependent of the beam waist radius wq 
(see Rcf. 15 and references therein). This beam repre- 
sents the light at the output of a fiber, its diameter 2wo 
being the diameter of the fiber core (about 8-10 /zm). 

Let us deal with the design of the spot-size converter. 
We assume that a TM polarized Gaussian beam (cen- 
tered at the origin of our coordinate system) propagates 
along the x-axis. This beam impinges a crystal consisting 
of a hexagonal array of Si cylinders in silica, the 2-axis 
being parallel to the cylinder axis. The size of this crys- 
tal along the y— direction has to be chosen slightly larger 
than the incident-beam-width to avoid flux escaping at 
the lateral borders of the crystal, hence we have chosen 
a crystal made of 26 rows of cylinders whose total length 
is about 3wo along the y— axis. Regarding its thickness 
along the x-axis, we used 13 columns of cylinders. The 
total number of cylinders is 318, which is conditioned 
by our calculation resources. This structure constitutes 
the lens-material (LM) of our integrated device. We con- 
sider a lattice constant ulm and radius r =0.294 ulm- A 
band structure calculation of the corresponding infinite 
crystal by means of a plane wave expansion method pre- 
dicts a forbidden gap for the TM-like modes in the range 
of frequencies 0.222-0.292 (in units of 2'Kc/aLM)- In or- 
der to have a low reflectance device we choose a working 
frequency of 0.197; i.e. below the first gap, where prop- 
agation is allowed in all directions. On the other hand, 
the condition imposed to the GA in order to produce 
the spot-size conversion is to maximize the x-component 
of the Poynting vector at one point (the focus) be- 
hind the cylinders. A full account of how the MST is 
used to get this magnitude can be found in Rcf. [15] and 
references therein. The CPU time needed to get the fit- 
ness of one individual is of the order of 6 s in a 2.8 GHz 
Pentium IV workstation. The coordinates of the focal 
point (xf,yf) are freely determined by us. Particularly, 
we impose the focal point to be located on the symme- 
try axis of the system (yj =0), and with a x— coordinate 
x/=4.60 Ao, where Aq is fixed by the working frequency 
(Ao = ai,M/0.197). Our results indicate that under this 



simple condition the losses due to reflectance diminish 
and, at the same time, it considerably reduces the waist 
of the propagated beam at the focus. 

We realize that the global maximum must correspond 
to a symmetric structure, which is obvious from the sym- 
metry of the problem. Then, the GA is forced to deal 
only with symmetric configurations. Therefore, the space 
search size is reduced from 2^ to 2^^+")/^ possible con- 
figurations, with N=318 and n=6, where N is the number 
of cylinders of the LM, and n the number of cylinders 
centered on the symmetry axis. In other words, from 
5.3x10^^ possible configurations to 5.8x10^^. We have 
employed a crossover operator with a probability C=0.4, 
which is maintained in the rest of calculations, and a pop- 
ulation P=400 individuals. The left panel in Fig. 1 shows 
the resulting spot-size converter device and the pattern of 
the corresponding electric field modulus. Two different 
scales are used to measure distances, lattice parameter 
units (bottom and left) and wavelength units (top and 
right). The spot-size conversion can clearly be observed 
in the right panel of Fig. 1, where Sx is represented along 
the direction parallel to the y— axis at the x— coordinate 
of the focal point (thin line), and for the incident beam 
at the y— axis (x =0) (thick line) . A fit of (x/, y) to a 
Gaussian along the y-direction gives a beam waist radius 
of Wo=0.30 Aq. This means that the spot-size converter 
has a conversion ratio of 3.30 Aq: 0.30 Aq « 11:1, suffi- 
cient for classical-PC waveguide mode coupling require- 
ments. Besides, this lens-type structure shows another 
property that is difficult to achieve by conventional lens 
design, its low F-number. According to Gaussian optics 
F=(27r)wQ/4Ao, which in our lens structure takes a value 
of only 0.47. On the other hand, the power loss, which 
is defined as the ratio of the transmitted power calcu- 
lated inside the focus waist (2wq) to that of the incident 
beam, takes the value of 1.0 dB. In other terms, a 79% 
of the incident power is squeezed and passes through the 
focus. At this point let us remain that Kosaka et al. [7] 
reported a photonic-crystal spot-size converter that re- 
duces the spot-size in the ratio 10:1, but no information 
is provided about the power passing trough the focus. 

Now, let us deal with the problem of designing an inte- 
grated device involving the coupling of the light squeezed 
by the spot-size converted into the PC waveguide. The 
parameters of the PC of what the waveguide is created, 
which we name guide-material (GM), must be chosen 
carefully in order to allow the coupling with a state local- 
ized in the waveguide. Firstly, the photonic band struc- 
ture has to have a full gap at the working frequency of 
the LM and secondly, a guided mode of the crystal with 
a missing row of cylinders, must exist at such frequency. 
Thus, we choose as the lattice parameter of the GM, 
o-GM=^-^0'LM- On the other hand, we keep the same ori- 
entation and cylinder radius as in the LM, and the PC 
waveguide is created along the FK-direction in the hexag- 
onal lattice. Figure 2 represents the TM-bands for the 
LM along FK as well as the projected band structure of 
the GM along the missing row of cylinders. It is shown 
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how the working wavelength Aq of the lens matches a 
mode in the guided band. 

If one simply places the entry of the waveguide at the 
focal point of the lens, one obtains a total insertion loss 
as high as 7.03 dB. This efficiency is calculated as the 
ratio of the total power that is transmitted through the 
waveguide to the incident power. The optimization pro- 
cess begins at this point by defining a set of 52 cylinders, 
which arc placed in front of the entrance in order to fa- 
cilitate the coupling to the waveguide mode, and letting 
the GA operate over it. These cylinders have the same 
symmetry, lattice parameter and orientation as the GM. 
The fitness was set equal to the sum of Sx calculated at 
30 points located along a transversal segment defined at 
the end of the waveguide; i.e., on the segment [-5.2 ulm 
, +5.2 ulm]- Now, each individual has only 32 genes due 
to symmetry reduction, and just 2^^ ~ 4.3x10^ possible 
combinations are available. The line (1) in Fig. 3 rep- 
resents the fitness of the best individual as a function of 
the number of generations. Each generation contains a 
population P of 400 individuals. The structure with the 
best fitness was achieved after 34 generations, the dashed 
horizontal line defines the maximum value obtained. A 
resulting power loss of 1.07 dB is obtained, which repre- 
sents a substantial advance in comparison with the total 
loss when no mouth-structure was considered. 

In order to reduce the insertion loss further, we allow 
the GA to act over the full structure made of the lens 
and the mouth. Now, its individual contains 162+32 
genes and, consequently, the space search size has been 
increased up to 2^^^+''^ ~2.5xl0'''**. In this case we con- 
sider that each generation contains a population of 600 
individuals. Figure 3 shows the result of running the 
GA on this structure, and Fig. 4 plots the optimized in- 
tegrated device (the wave-guide coupler) together with 
its electric field modulus pattern. The insertion loss 
predicted for this structure is as low as 0.61 dB. This 
means that 87% of the impinging light passes through 
the waveguide and is detected at the output. In fact, 
this value is underestimated since it does not include the 
light reflected at the end of the waveguide by finite size 
effect. Therefore, the coupling efficiency predicted by 
this new structure is comparable with the ones reported 
in the literature. Thus, a waveguide-to-fiber coupling im- 
provement exceeding 2dB per converter is shown in Ref. 
[2]. References [3] and [4] reported two different tapered 
couplers that numerically have over 90% power transmis- 
sion. The simulations by P. Sanchis et al. [5] reported 
that a transmission over 84% can be achieved if defects 
are put on a planar photonic crystal tapered waveguide. 
Finally, the J-coupler proposed by Prather et al. [6] pre- 
dicts a coupling efficiency of 91%. In spite of the large 
coupling efficiencies reported for these devices, in com- 
parison with the one here designed, the low compactness 
and integrability of these devices are their main draw- 
backs when they are compared with the one designed by 
evolutionary programming. 

Although our simulations involved the simplifying as- 



sumption of 2D PC's, it should be noticed that such 2D- 
periodic crystals can be studied in actual 3D crystals. 
For example, our recent simulations of 2D prism-^^ re- 
produces fairly well the behavior of the actual structures 
in the microwave regime. We conjecture that similar re- 
sults might also be obtainable in the optical regime by 
using a PC-slab sandwiched between multilayers films 
with a large gap. Also, the preceding discussion focused 
on the TM modes of a structure based on "dielectric- 
scatterers-on-background" . However, based on the gen- 
eral method presented here similar devices based on 
"holes- in-dielectric" structures can also be obtained. 

In summary, this work has shown that a genetic al- 
gorithm used in combination with multiple scattering is 
able to design photonic crystal-based structures with spe- 
cific optical properties. By using this approach, we have 
demonstrated the viability of fabrication of waveguide 
couplers with low losses and a high spot-size conversion 
ratio. Besides, they can be integrated in the same wafer 
with other planar structures in a simple single-step litho- 
graphic process. It can be said that this method repre- 
sents a substantial advance in order to reach the ultimate 
goal in optical devices design, that is molding the flow of 
light. 
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FIG. 1. (Left panel) Focusing effect produced by a spot-size 
converter (Ions) based on a 2D photonic crystal (black dots). 
The pattern of the electric field modulus is represented in 
a wide spatial region. The lenght scales are given in terms 
of the lattice parameter employed in the design of the lens, 
ulm, as wells as in terms of the working wavelength of the lens 
Ao- (Right panel) The x-componcnt of the Poynting vector 
represented at the focal point (a:/=4.60 Ao). A scale color 
is used, red (blue) color means maximum (minimum) electric 
field modulus. 



FIG. 2. (Left panel) Photonic band structure of the pho- 
tonic crystal employed as lens material (LM) along the FK 
direction (see inset). (Right panel) Projected band structure 
of TM modes in the guide material (GM). The grey regions 
represents the continuum of extended states in the photonic 
crystal. The white regions are the photonic band gaps. The 
waveguide is formed by removing one row of Si rods as shown 
in the inset. The black line defines the guided modes inside 
the waveguide. The frequencies in both panels arc given in 
reduced units of the respective lattice parameters, olm and 
acM- The horizontal dashed line defines the working fre- 
quency of optical devices under design. 



FIG. 3. Running the genetic algorithm 
on the lens+waveguide integrated structure. Line (1) defines 
the fitness of the best structure when the optimization only 
acts over the cylinders on the mouth of the waveguide, the 
dashed line defines the maximum fitness achieved. Line (2) 
represents the corresponding result when the GA acts over 
the cylinders in the lens and mouth simultaneously. Line (3) 
being the averaged fitness over the total population in this 
case. The vertical lines separate the regions where different 
values of mutation parameter are employed. 



FIG. 4. The optimized wavcguided-coupler device (dots) 
obtained by a genetic algorithm. The orange rectangles en- 
closed the cylinders where the GA is applied. A scale color is 
used; red (blue) color defines the maximum (minimum) elec- 
tric field modulus. 
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